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Abstract 


The  object  of  this  research  was  to  compile  spectroscopic  data  for  the 
strongest  bands  of  the  molecules  NO,  NH^  and  HNO^  and  to  put  it  in  the  format 
of  the  AFGL  Atmospheric  Absorption  Line  Parameters  compilation.  This  has 
been  done  for  the  1-0,  2-0,  and  3-0  bands  of  NO  (molecule  8).  For  the  inversion 
spectrum  of  NH^  (molecule  11);  for  the  pure  rotation  spectrum,  and  the  bands  at 
7.5y  and  11. 3u  of  nitric  acid  (molecule  12).  These  are  all  trace  gases  in  the 
Earth's  atmosphere.  For  NO  and  NH^  papers  have  been  submitted  for  publication 
(see  appendices)  while  for  HNO^  a description  of  the  data  are  of  sufficient 
accuracy  to  permit  a positive  identification  of  the  trace  gas  (pollutant) 
molecule. 


ACCESSION  ttr 


NTiS 

DOC 

VNANNOinCED 

JUSTIFICATION.. 


Wlilte  Sactia 
l«tf  SKilM  □ 
□ 


IT 


DISTRIBUTION/AVAlUaiun  COOES 


DM. 


AVAIL  Miil/or  SPECIAL 


A 


I \ 


i 1 


An 


Improved  Fit  to  the  Inversion  Spectrum  of 


L.  D.  G.  Young  and  A.  T.  Young 
Physics  Department 
Texas  A&M  University 
College  Station,  Texas  77843 


Abstract 

A Fade  approximation  with  19  terms  is  used  to  fit  to  measured  line  positions 
in  the  inversion  spectrum  of  ammonia.  It  gives  an  order  of  magnitude  better 
fit  than  an  exponential  (Costain)  function  with  21  terms,  and  two  orders  of 
magnitude  better  fit  than  a polynomial  with  21  terms. 

Introduction 

In  the  course  of  preparing  data  for  the  AFGL  trace  atmospheric  gas  line 
parameter  compilation,  we  discovered  that  Taylor's  published^^^  line  positions 
for  the  vg  fundamental  band  of  NH^  at  10. 5u  could  be  off  as  much  as  -1.5  cm"\ 
or  15,000  times  the  estimated  error  in  measured  line  positions^^’^^ . As  a 
result,  measured  line  pocitions  were  used  in  the  data  compilation^^*^^ , and 
we  began  to  look  into  the  reason  for  the  discrepancy  between  the  measurements 
and  the  calculations. 

One  of  the  first  things  we  did  was  to  loOk  at  the  accurate  microwave 
measurements  of  Poynter  and  Kakar^®^.  They  could  not  fit  their  own  measurements 
of  the  inversion  spectrum  which  "are  good  to  ±0.005  MHz,”  with  a standard 
deviation  less  than  0.034  MHz.  They  tried  two  approaches:  First  they  assumed^^^ 
that  their  data  could  be  fit  by  power  series  in  J(J  + 1),  the  square  of  the  total 
angular  momentum,  and  K,  the  projection  of  the  angular  momentum  on  the  symmetry  axis 


of  the  molecule.  That  is,  the  inversion  energy  could  be  represented  by 


E = BJ(J+1)  + (C-B)K^  + Dj  J^(J+1)^  + J(J+1)K^ 

+ + Hj^^J(J+l) 

* Hjjj  jVi)3  * K«  . . . . (1) 

P 

with  the  line  positions  being  given  by  a similar  power  series  in  J(J+1)  and  K . 
For  a polynomial  of  this  type,  Poynter  and  Kakar  found  10  terms  gave  a standard 
deviation  of  a » 23.14  MHz,  15  termsgaveo  = 2.308  MHz,  and  21  terms  gave 
o = 0.228  MHz. 

Next  Poynter  and  Kakar  tried  the  approach  suggested  by  Costain^®^  in  a 
paper  entitled  "An  empirical  formula  for  the  microwave  spectrum  of  ammonia". 
Here  the  inversion  frequency  is  represented  by 

Y = a exp[bJ(J+l)  + cK^  + dJ^(J+l)^  + eJ(J+l)K^  + fK^  + . . .]  (2) 

Using  a power  law  within  the  exponent  gave  Poynter  and  Kakar  a somewhat  better 

fit  for  10  terms;  a = 0.495  MHz;  however  they  could  only  get  an  order  of 

magnitude  improvement  in  going  to  21  terms:  o = 0.039  MHz. 

Clearly  neither  of  these  approximations  could  approach  their  experimental 

accuracy  without  a major  increase  in  the  number  of  constants  being  fit. 

(91 

In  referring  to  work  on  the  NO  molecule,  Mizushima'  ' noted  that  high 
order  perturbations  can  lead  to  matrices  in  the  form  of  a Pade  approximant. 

Thus  we  decided  this  might  be  a useful  approach  for  an  empirical  fit  to  the 
inversion  spectrum  of  NH^. 

We  first  tried  fitting  the  inversion  energy  with  a fraction  in  which 
both  numerator  and  denominator  consisted  of  polynomials,  such  as  given  in 


eqn.  1.  While  this  approach  gave  a better  fit,  for  the  same  number  of 
constants,  than  either  of  the  expressions  for  the  inversion  energy  which 
Poynter  and  Kakar  had  used,  we  found  we  could  get  an  even  better  fit  if  we 
used  the  following  expression 


E « [a  + bjF(J,K)  + b^K^  + djjF(J,K)^  + dj^F(J,K)K^  + d^^K^ 
[1  + b*F(J,K)  + b*K^  + dj^F(J,K)K^  + d*,^K^  + djjF(J,K)^ 


where  F(J,K)  = [r  - J(J+1)]. 

Lines  with  K *=  3,  6,  9,  ...  need  to  have  an  additional  J-dependent  term 
added,  since  they  are  affected  by  a K-type  splitting^^®’^^^.  Thus  for  K = 3 
one  must  add 


AE3  = (-1)'^  f3(J)[ai  + a2J(J+1)]/[l  + a3J(J+l)] 


where  f3(0)  = J(J+1 )[J(J+1 ) - 2][J(J+1)  - 6].  For  K = 6,  the  correction  is' 


AEg  = fg(J)C6i  + 02'^(J+1)] 


where  fg(J)  « f3(J)[J(J+l)  - 12][J(J+1)  - 20][J(J+1)  - 30].  For  K - 9,  the 
correction  is 

AEg  - (-1)'^  fg{j)Y^ 

where  fg{J)  * fg(J)[0(J+l)  - 42][J(J+1)  - 56][J(J+1)  - 72].  Using  the 


constants  given  in  Tables  1 and  2,  we  found  we  could  fit  Poynter  and  Kakar's 
measurements  with  a standard  deviation,  for  a line  of  unit  weight,  of  6.9  kHz, 
which  is  close  to  their  nominal  accuracy  of  5 kHz.  Figure  1 illustrates  how 
the  standard  deviation  drops  as  the  number  of  parameters  increases. 

Poynter  and  Kakar's  poor  fit  was  due  to  three  factors:  First,  they  in- 
cluded several  lines  of  lower  accuracy  from  the  NBS  compilation^^^^ . Second, 
they  gave  all  lines  equal  weight,  even  though  they  knew  some  lines  were  less 
accurate.  When  observations  of  unequal  quality  are  combined,  they  should  be 
given  weights  inversely  proportional  to  the  square  of  the  error  for  each  set 
of  data.  Third,  they  did  not  reject  any  bad,  or  discrepant,  data. 

Experience  has  shown  that  any  iterative  least-squares  program  must  be 
able  to  reject  aberrant  data.  As  a practical  matter,  we  have  found  that 
good  results  can  be  obtained  if  we  reject  data  having  a residual  greater  than 
4a  (four  standard  deviations).  We  found  the  (12,1),  (11,8)  and  (6,2)  lines 
to  have  residuals  of  -339,  -108  and  -37  kHz,  respectively;  these  are  re- 
jected by  our  4-a  criterion.  ' 

» 

Although  the  observed  line  positions  are  fit  to  a few  kHz,  the  least-  ! 

squares  values  of  the  parameters  in  Tables  1 and  2 have  rather  large  errors.  j 

This  is  primarily  due  to  the  large  correlation  between  numerator  and  denominator  i 

* ( 

terms.  Thus,  it  is  necessai^y  to  give  the  coefficients  several  figures  beyond  t, 

the  leading  significant  digits  in  their  standard  errors.  One  readily  sees  J j 

that,  to  retain  1 kHz  accuracy  up  to  quantum  numbers  on  the  order  of  15,  we  { { 

must  specify  b's  to  about  10~^,  d's  to  about  10"^,  and  h’s  to  about  10"^®.  i 

We  have  been  careful  to  give  enough  significant  figures  to  allow  others  to 
reproduce  our  results. 

The  improved  accuracy  of  our  fit  has  allowed  us  to  more  accurately  model 
the  perturbations  which  affect  the  K = 3 and  K = 6 lines,  and  to  detect,  for 
the  first  time,  the  perturbation  of  K = 9 at  a signal/noise  ratio  of  3.4 


standard  deviations. 


Table  give  the  computed  frequencies  for  the  ammonia  inversion  lines  in  the 


microwave  spectrum. 
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Table  1.  Coefficients  in  the  Pade  (rational)  approximation  to  the 
inversion  spectrum  of  amnwnia. 


Numerator  Coefficients 

(HHT) 


Denomi nator  Coefficients 
rSTniehsionless) 


23785.8985  ± 0.0038 


'J 

79.750378i3.3871 

* 4 

bjXlO^ 

-30.181  1.42 

'k 

21.145474±1.7309 

bj^xlO^ 

-16.31  1 0.73 

10.355615±0.9322 

2.32810.516 

Ij^xlO^ 

6.524901^0.4487 

d*,xlO« 

4.59810.210 

l^^xlO 

- 0.50860810.2535 

d*,xl0‘ 

1.11310.077 

54.215901  7.80 

* 10 

5.561  5.87 

58.94479  1 6.85 

* 10 

-45.931  3.48 

- 4.744461-1.45 

♦ in 

hji^l^xlO 

-21.451  1.34 

- 8.635031  0.97 

- 3. 58 1 0.33 

36124.546  ± 3.64 

MHz 

206.4072 ± 18.3 

MHz 

77383.461  ± 645. 

(dimensionless) 

16.562  ± 1.60 

. MHz 

17.896  ± 7.76 

MHz 

15.76  ± 4.61 

MHz 

Table  3.  Computed  inversion  frequencies,  in  MHz. 
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Figure  1 


. Standard  deviation  for  the  fit  to  the  inversion  spectrum  of  ammonia 
as  a function  of  the  number  of  parameters  used  in  the  fit. 


Nitric  Acid  Vapor  (Molecule  12) 


The  nitric  acid  molecule  is  approximately  an  oblate  symmetric  top,  since 
the  ground  state  rotational  constants  areA^=0.4340  cm"\  8^=0.4036  cm'^  and 
Cq*0.2088  cm"^.  The  C axis  is  perpendicular  to  the  plane  of  the  molecule, 
while  the  dipole  moment  lies  in  the  plane.  Note  that  ABiBftJ2C.  This  gives 
rise  to  a rotational  spectrum  which  is  interesting,  since  a large  number  of 
coincidences  of  line  positions  occur.  It  results  in  a very  simple  spectral 
pattern  under  moderate  resolution.  One  can  write  the  energy,  E,  as 

E * 0.41865  J(0+1)  -0.20985K^  + Smaller  terms 

(Flemming,  1974).  This  spectrum  has  the  general  appearance  of  that  for  a 
diatomic  molecule  with  a line  spacing  of  -v  0.4  cm~\ 

The  line  positions  are  given  by  Flemming  and  Wayne  (1975).  They  published 
a spectrum  taken  between  10  and  40  cm”^  with  a resolution  of  0.05  cm"^.  The 
pressure  was  4.5  torr  and  the  pathlength  was  1 meter.  Using  the  values  of  the 
halfwidth  for  self  broadening  (o^  = 0.8  cm”^)  and  nitrogen  broadening  (a^=0.13  cm”^) 
as  reported  by  Brockman  et  a1.  (1978),  and  assuming  a square  root  absorption  law, 
we  have  derived  values  for  the  line  intensity.  The  portion  of  the  7.5w  band 

measured  in  the  laboratory  by  Fontanella  et  a1.  (1975)  for  a total  pressure  of 

17  18  3 

66  mb  and  amounts  of  HNO^  vapor  ranging  from  0.27  x 10  to  0.44  x 10  molecules/cm  , 

has  also  been  used  to  derive  line  intensities  for  that  band.  It  also  appears  to 

have  a "line"  spacing  of  ''-0.4  cm”^.  The  lower  resolution  data  of  Goldman  et  al. 

(1971)  and  Moskalenko  (1972)  have  not  been  used,  as  the  individual  spectral  lines 

are  not  resolved.  However,  the  integrated  intensity  of  the  7.5u  band,  obtained  by 

Goldman  et  al.,  was  used  to  check  the  line  intensities  derived  from  the  spectra 

published  by  Fontanella  et  al.  Approximate  values  of  the  rotational  quantum  numbers 

are  given  in  the  AFGL  trace-gas  compilation.  The  line  strengths  for  both  the  pure 

rotational  spectrum  and  the  7.5m  band  are  estimated  to  be  accurate  to  i 20%. 


Brockman  et  al.  (1978)  have  made  very  high  resolution  measurements  of  a 
portion  of  the  11. 3p  band.  No  quantum  numbers  have  been  assigned  to  the  1080 
lines  measured.  They  estimated  the  uncertainties  of  the  line  strengths  as  ± 20%, 
and  of  the  line  positions  as  ± 0.005  cm"\  The  pure  rotational  line  positions 
have  an  estimated  precision  of  ± 0.003  cm"\  while  the  positions  of  the  lines  In  the 
7.5u  band,  are  considerably  less  certain,  having  been  read  off  of  the  published 
spectra;  they  are  only  accurate  to  i 0.1  cm“^. 

In  putting  the  band  Identification  on  the  IBM  cards,  we  ran  Into  the  problem 
of  not  enough  room  in  the  card  format  for  all  the  quantum  numbers  to  be  listed  In 
the  usual  way.  I have  given  the  band  location  In  microns  and  the  lower  (ground) 
state  simply  represented  as  0. 
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Abstract 

The  calculation  of  line  positions,  energy  levels  and  line  intensities 
for  the  NO  fundamental,  done  for  the  AF6L  atmospheric  trace  gas  data  tape, 
is  presented.  A simultaneous,  least-squares  fit  to  measured  frequencies  for 
the  fundamental  and  satellite  bands  of  to  its  microwave  spectrum,  and 

to  two-photon  spectra  in  the  microwave  region  has  yielded  new  spectroscopic 
constants.  The  eigenvectors  obtained  were  used  to  compute  line  intensities 
for  these  bands.  Line  halfwidths,  based  on  a linear  fit  to  the  measurements 
of  Abels  and  Shaw,  were  used  in  the  data  compilation. 


^Participant  USAF-ASEE  Summer  Faculty  Research  Program,  at  the  Air  Force 
Geophysics  Laboratory,  1976. 
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Introduction 


There  has  been  a group  of  papers  recently  on  the  NO  fundamental . 

This  interest  derives  from  the  facts  that  nitric  oxide  is  a trace 

(pollutant)  gas  in  the  earth’s  atmosphere; and  is  also  one  of  a few 
2 

molecules  with  a n electronic  ground  state  and  hence  is  inherently  in- 
teresting to  spectroscopists. Our  calculations  differ  from  other 
recently  published  results  since  we  have  included  microwave  data  in  addi- 
tion to  data  on  the  fundamental.  We  have  also  included  higher-order  terms 


(e.g.  H^)  in  the  expression  for  the  rotational  energy  levels. 


Recent  high  resolution  measurements  by  Valentine  et  al. 


(10) 


on  the  in- 


frared spectrum  of  NO  have  yielded  substantially  more  accurate  line  positions 


for  the  two  fundamental  bands  (^n 


1/2 


satellite  bands  (^n 


1/2  ^ "3/2  ^"3/2 


> 2 2 

'11^2  ^3/2  ^3/2^  their 

■n^yg)-  addition  to  these 


measurements.  Pine  and  Nill  have  recently  measured  (but  not  yet  published) 
lambda-doubling  in  the  fundamentals,  using  two-photon  spectroscopy  which 
has  the  accuracy  inherent  to  microwave  spectroscopy.  By  combining  these  new 
data  with  the  older  microwave  measurements on  the  ground  states  [the 
l/2->3/2  and  3/2->5/2  lines  measured  by  Gallagher  and  Johnson^^^^  and 
the  ^^2/2  ~ measured  by  Favero  et  al.^^^^  ],  we  believed  that 

accurate  eigenvectors  and  energy  eigenvalues  could  be  obtained  for  both  the 
V = 0 and  the  v = 1 states.  Our  data  base  consisted  of  298  measured  fre- 
quencies, which  were  weighted  according  to  the  accuracy  of  the  experi- 
mental measurements.  The  weights  used  in  the  computer  program  are  based 
on  the  reciprocal  of  the  squares  of  the  stated  errors;  lines  having  residuals 


greater  than  four  standard  deviations  were  flagged  and  rejected  for  the 

next  iteration.  In  our  original  calculations,  we  included  the  lines  in 

the  pure  rotational  spectrum^ and  older  measurements  of  the  funda- 

(20'-29) 

mental  band  in  the  infrared  • 

Line-Frequency  Calculations 

2 

The  appropriate  theory  of  n electronic  ground  states  dates  back  to 

(301  22 

Van  Vleck.  ' ' For  NO,  the  two  ground  states,  11^2  ^‘i/2 

turbed  by  an  upper  electronic  state.  Keck^^^^  and  James^^®^  have 

indicated  the  matrix  elements  for  the  appropriate  two  3x3  matrices  of 
* 

the  Hamiltonian,  xiz. 


uS,a 

Hi  2 

”l3 

- 

H^  + A (L-S) 

= [ 

H22 

H23  ] 

(1)  ’ 

”31 

«32 

^33 

1 

• 1 

Here  the  superscripts  £ and  a^  refer  to  the  symmetric  and  antisymmetric 
basis  functions;  also. 


I 


”li*  * 1 (J  + 1/2), 

(2a) 

’ r 

I 

1 

Hfl®  = a^  + 2bg  + 2bg  (0  + 1/2), 

(2b) 

■ i 

t 

1 

H^3  = 2bg  [ (0  + 3/2)(J  + 1/2)  ] ^ 

(2c) 

1 

*. 

“22  “ % ■ * V2)^(B^  - Aj)  - (J+l/2)\  + (J+l/2)\  , 

(2d) 

H23  - [B^  + 2D^  - 2Dv(J  + 1/2)^  3 [(J  + 3/2)  (J  - 1/2)  ] , 

(2e) 

“33  “ % * - 2(Bv  + 20J  + (J  + 1/2)^  (B^  + 40^  + Aj)  - 

(0  + l/2)\  + (J  + l/2)\  , 

(2f) 

it  ~ - - - ^23) 

See  also  Dousmanis,  Sanders  and  Townes.'  '0 

j 
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where  the  upper  sign  in  egn  (2d)  and  (2b)  refers  to  the  £ matrix  element  and 
the  lower  sign  refers  to  the  a matrix  element.  Here  the  constants  G , 

— V 

B^,  D^,  and  have  their  usual  meaning,  Aj  allows  for  the  J dependence 
of  the  spin-orbit  interaction  constant  A,  and  A^  allows  for  its  depen- 
dence on  vibrational  state.  The  A-doubling  parameters  are  a^  = <111 A LylE> 
and  bg  = - <7rlB 

Holding  the  constants  and  Bj.  fixed,  and  using  the  spectroscopic 

constants  for  the  fundamental  reported  by  Olman  et  al.^^^^  as  starting 

values,  we  proceeded  to  do  an  interative  least-squares  correction  of  the 

parameters  given  in  Table  1.  The  measured  line  frequencies  and  doublet 

splittings  were  compared  with  frequencies  computed  from  the  difference  of 

energy  eigenvalues  between  the  upper  and  lower  states.  The  constants  a^ 

and  bg  were  poorly  determined  if  the  microwave  measurements  of  lambda 

doubling  in  the  fundamentals  were  not  included.  Giving  these  measurements 

a weight  of  500  (compared  to  a weight  of  1 for  the  infrared  data)  resulted 

in  well -determined  values  for  a„  and  b . The  values  of  and  „ could 

e e E z-n 

2 

be  varied  by  five  percent  with  a negligible  effect  on  the  values  of  the  n 
state  spectroscopic  constants  obtained  on  fitting  the  data.  The  standard 
deviation  found  for  the  frequency/lambda-doubling  fit  was  0.0011  cm~^  for 
a measurement  of  unit  weight.^ 

Originally, the  measurements  of  Keck^^^^  and  James^^®^  were  included. 

Both  of  these  sets  of  data  had  substantially  higher  residuals  than  the  data 
of  Valentine  et  al.^^®^  It  was  decided  that  James*  and  Keck*s  data  were 
contributing  mostly  "noise"  to  the  solutions,  so  these  data  sets  were  dropped. 

^The  microwave  lines  have  a typical  weight  of  7000;  Pine's  measurements  of 
lambda  doubling  have  weight - 1000;  Valentine  et  al . measurements  have  weight 
2 except  for  "bad  lines"  which  have  weight  0.5. 
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Recent  measurements  of  19  NO  lines  made  using  a CO  laser  given 

weights  comparable  to  those  of  the  273  lines  measured  in  the  infrared.  The 

pure  rotation  lines  measured  by  Palik  and  Rao  and  by  Hall  and  Oowling^^^^ 

had  large  systematic  deviations  and  were  dropped  from  the  final  solution.  The 

relatively  large  errors  and  hence  low  weights  for  these  lines  leads  us  to 

n 5l 

believe  that  Rabach'e  calculations  ' ' are  inaccurate,  since  he  used  the  pure 

rotational  lines  in  his  fit  to  the  molecular  constants  for  NO.  In  addition, 
his  calculations  assume  the  zero  energy  level  lies  midway  between  the 
^112^2  states  and,  therefore,  his  computed  energy  levels  are  inaccurate  by  at 
least  '\/2. 

Table  1 gives  the  values  obtained  for  the  spectroscopic  constants.  The 

constants  of  Valentine  et  al.,^^^^  using  the  theory  of  Albriton  and  Zare,  are 

given  for  comparison.  With  the  exception  of  H^,  we  have  a smaller  standard 

deviation  for  our  constants  than  the  latter  have.  Finally,  it  should  be  noted 

that  the  use  of  effective  constants^^^^  in  fitting  Keck's  measurements of 

the  fundamental  yield  values  with  considerably  larger  errors,  but  taking  that 

* 

into  account,  give  results  in  agreement  with  the  present  calculations. 

By  Including  the  very  accurate  microwave  data,  we  have  obtained  much 
greater  accuracy  for  the  low-order  rotational  constants  and  the  A-doubling 
constants. 


Line  Intensity  Calculations 

When  we  diagonalize  the  Hamiltonian  matrices,  we  obtain  the  eigenvec- 
tors V,  J),  in  addi- 

2 

tion  to  the  energy  eigenvalues.  Using  the  n eigenvectors,  it  is  straight- 
forward to  compute  the  line  strengths  using  the  electric  dipole  moment 

(331 

transition  probabilities  for  symmetric  tops.'  ' For  the  two  fundamental 


The  results  are  indicated  in  Table  2. 


~~T^ 


\ 


I 


4 


t i 
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2 2 2 2 

bands,  11^2  “ ^^2  ^3/2  " 1^3/2  • obtain  essentially  the  same  line 

(34) 

strengths  as  James.  We  compared  our  calculated  line  strengths  with  James' 
for  the  two  satellite  bands.  For  the  satellite  bands,  James  only  published 
line  strengths  for  the  ^2/2  “ ^1/2’  ^'^^osition.  Our  present  calculations 

yield  somewhat  different  values  for  the  line  strengths.  We  can  reproduce  James' 
results  by  setting  the  lambda-doubling  constants  a^  and  b^  equal  to  zero  and 
using  the  spectroscopic  constants  we  obtained  for  the  ground  state  as  the  con- 
stants for  both  the  upper  and  lower  vibrational  states.  Goldman  and  Schmidt^®^ 
noticed  a similar  result  in  their  calculations  as  compared  with  James. 

Line  Halfwidths  Used  in  the  AFGL  Data  Compilation 

(351 

A comparison  of  line  widths  measured  by  Hochard-Oemolliere  et  al.'  'by 
Abies  and  Shaw^^^^  with  the  linewidths  calculated  by  Tejwani  et  al.^^^^  show 
considerable  discrepancies.  In  our  judgement,  a linear  fit  to  the  measurements 
of  Abies  and  Shaw  should  give  reasonably  good  values  for  atmospheric  trans- 
mission calculations. 

Discussion 

By  diagonalizing  the  Hamiltonian  given  in  eqn.  (1)  and  fitting  the 

computed  energy  eigenvalues  to  experimentally  observed  line  positions,  we 

have  obtained  quite  precise  energy  values  of  the  v"  = 0 and  the  v'  = 1 states 

of  In  addition,  by  using  the  eigenvectors  obtained  in  the  above 

process,  we  have  computed  accurate  line  strengths.  Valentine  et  al.^^®^ 

(141 

and  Johns  et  al.'  ' have  made  calculations  very  similar  to  ours.  We  give 
a comparison  of  the  results  in  Table  1.  This  shows  that  our  values  for  the 
constants,  as  a whole,  have  values  intermediate  to  those  found  by  these  other 
authors.  The  values  of  a^  and  b^  were  computed  from  their  published  values 
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of  p and  q and  E as  given  by  Herzberg^^^^  The  discrepancy  between  our  calcu- 

I*  11 

(13) 

lations  and  those  of  Valentin  et  al.'  'for  values  of  H^and  Aj  is  puzzling, 
since  these  constants  are  determined  from  the  high  J lines  and  we  have  employed 
the  same  data.  The  differences  in  the  other  constants  can  be  attributed  to 
our  inclusion  of  measurements  in  the  microwave  region. 

The  calculations  of  Meerts^^^^  and  Kristiansen^^®^  employed  a different 
perturbation  scheme  and  a different  set  of  spectroscopic  parameters.  The 
simplest  way  to  compare  these  results  with  ours  is  in  terms  of  effective 
constants.  This  is  done  in  Table  2.  Kristiansen^^®^  used  the  data  of  Keck^^^^ 
and  of  James  and  Thibault^^®^  for  infrared  data;  he  used  those  of  Hall  and 
Dowling^^^^  for  the  rotational  spectrum,  and  the  same  data  as  we  used  for  the  micro- 
wave  lines.  In  Kristiansen's  determination  of  line  parameters  for  NO,  the 
spectroscopic  constants  0^  and  were  not  included  in  his  fit  of  16  parameters 
to  159  lines. 

Our  calculations  included  these  parameters,  and  we  fit  to  298  lines 
and  lambda  splittings.  Meerts^^®*^®^  fit  the  hyperfine  splittings  of  NO  lines 
in  the  microwave  region,  using  91  observed  transitions  in  the  frequency  range 
from  0.6  MHz  to  258  GHz.  He  required  22  molecular  constants  to  describe 
the  observations. 

While  our  constants  a^  and  b^  show  a slight  dependence  on  the  vibrational 
state,  one  would  not  have  expected  this.  However, when  these  constants  are  set 
equal  to  the  mean  value  for  the  two  vibrational  states,  the  rms  error  in  a 
line  position  of  unit  weight  increases  from  a * 0.000647  to  0.005241.  In 
addition, the  program  rejects  all  of  Pine’s  measured  A splittings  in  the  R- 
branch  of  the  1/2-1/2  subband.  Thus, the  difference  between  the  a^  (and  b^) 
values  for  the  v = 0 and  v * 1 states  is  both  statistically  significant 
and  well  determined. 


The  input  data  for  the  AFGL  atmospheric  trace  gas  tape  should  be  accurate 
to  better  than  0.0004  cm”\  which  is  about  the  accuracy  routinely  available 
in  CO  laser -line  data. 
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Table  I.  Spectroscopic  constants  for  N 


present  results 
a =,000647 


Valentine 


Johns  et  al. 


(13)  (14) 


B.  1.69613891 
±.00000116 


1.6961360 
± .0000032 


1.6961087 
± .0000040 


D 5.49149  X 10" 
±.00374 


5.47498  X 10 
±.00055^' 


5.495  X 10. 
±.013 


H „ 2.072  X 10 
±.258 


A 122.94865 
.00145 


1.75852  X 10 
.00332 


90.130 

±.352 


bnn  -71564 
±.00279 


44.000 


Bj-o  1.9870 


G ,1875.97116 
±.00041 


B„, 1.67856953 
.00000121 


0 ,5.50779  X 10' 
.00333 


H .1.837  X 10” 
±.216 


A„, 122. 70577 
±.00304 


A., 1.70448  X 10' 
.00304 


a., 90. 236 
±.276 


1.346  X 10  “ 
±.300 


123.02774 

±.00011 


3.47621  X 10" 
.00050 


0.7130 

±.0004 


44.0003 


1875.97240 

±.00007 


1.6785673 

±.0000030 


5.49278  X 10 
.00050 


1.209  X 10" 
±.255 


122.78249 

±.00011 


3.36993  X 10' 
.00050 


122.8778 

±.0490 


-9.89  X 10' 


116.3b 
+ 0.3 


44.0003 


1875.9725 

±.0010 


1.6785438 

±.0000107 


5.513 

+.013  X 10" 


122.6427 

±.0440 


9.91  X 10' 


115.7 

±0.3 


These  values  were  converted  to  the  notation  used  in  the  present  paper. 
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v=2  and  v=3  States  of  Nitric  Oxide 
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Abstract 


The  calculation  of  line  positions,  energy  levels  and  line  intensities 
for  the  NO  overtone  bands  2-0  and  3-0,  done  for  the  AFGL  atmospheric  trace 
gas  data  tape  is  presented.  A simultaneous  least  squares  fit  has  been  made 
to  the  measured  frequencies  of  the  overtone  bands,  using  previously  determined 
constants  for  the  ground  (v=0)  state.  This  has  resulted  in  new  spectroscopic 
constants  for  the  v=2  and  v=3  states.  The  eigenvectors  were  used  to  compute 
line  intensities,  and  the  same  line  halfwidths  as  were  used  for  the  fundamentals 
are  used  for  the  overtone  bands. 


Introduction 


In  an  earlier  paper'*'  we  reported  the  results  obtained  from  a least- 
squares  fit  to  fundamental  band,  and  the  details  of  the  line  frequency 
fit  are  presented  there.  The  same  computer  program  was  used  for  the  over- 
tone bands.  Once  again,  we  fit  to  the  recent  laboratory  measurements  of  NO 
line  positions  given  by  Henry  et.a1.'  . These  latter  authors  used  the  method 

of  Zare  et.  al.^^^  to  obtain  spectroscopic  constants  for  the  three  bands  (1-0), 
(2-0)  and  (3-0).  However,  for  these  bands,  they  only  included  lambda  doubling 
of  the  lines  as  a correction  to  the  unsplit  line  positions.  Our  computer 
program  includes  A-doubling  explicitly. 

The  Overtone  Bands  for  the  (2-0)  band,  143  lines  have  been  published  by 
Henry  et  a1.  Blends  have  been  given  a weight  of  zero  in  our  calculations. 

All  other  lines  have  a weight  = 1.0.  Our  fit  to  their  measured  frequencies 

-4  -1 

gives  a standard  deviation  of  3.7x10  cm  , while  Henry  et. al.  report  a 

-4  -1 

standard  deviation  of  5x10  cm  . 

There  is  a reason  for  our  improved  fit  to  their  data:  Henry  et  al.  were 

forced  to  assume  that  the  lambda -doubling  parameters  p and  q (analgous  to 

our  a„  and  b^)  had  the  same  values  for  the  vibrational  states  v=l,  v=2,  and 
e e 

v»3.  We  don't  have  to  make  this  assumption  for  the  (2-0)  band. 

For  purposes  of  comparison  with  the  data  of  Henry  et  al.,  we  did  assume 

the  v*l  constants  could  be  used  for  the  (2-0)  band.  This  produced  a larger 

standard  deviation,  4x10"^  cm~^ , than  our  fit  where  both  a^  and  b^  were 
allowed  to  vary.  For  this  band,  both  a^  and  b^  were  fairly  well  determined. 

For  the  (3-0)  band,  Henry  et  aU  only  published  frequencies  for  112  lines 
and  lambda  splittings.  Here,  they  only  measured  two  lines  in  the  Q branch 
of  the  with  lambda  doubling.  This  simply  is  too  few  data  for 


our  program  to  find  both  and  b^.  Since  is  simply  a constant  in  our 
Hamiltonian,  we  can  hope  to  determine  it.  The  constant  b^  appears  as  a 
coefficient  of  terms  involving  the  rotational  quantum  number  J and  the  J 
values  of  the  split  lines,  for  the  band,  are  only  J = 3/2,  and 

J = 5/2.  As  a result  we  can't  get  a good  value  for  b^.  We  have  assumed 


values  for  b^  fixed  at  0.74,  0.75  and  0.76  cm' 


The  three  fits  to  the 


-4  -1 

frequencies  measured  all  had  a standard  deviation  of  a = 5.81x10  cm  , 

so  we  present  our  results  in  Table  1,  only  for  b = 0.75.  The  constants  of 

(21  ^ -4 

Henry  et,al.'  , also  shown  in  Table  1,  had  a standard  deviation  of  a = 7x10  cm. 

We  believe  our  smaller  standard  deviation  resulted  from  including  H^  terms  in 

the  rotational  energy  expression  and  our  use  of  smaller  values  for  the  lambda 

doubling  parameters. 


References 


L.  D.  G.  Young,  A.  J.  Young,  S.  A.  Clough,  F.  X.  Kneizys,  JQSRT  (in  press). 
A.  Henry,  M.  F.  LeMoal , Ph.  Cardinet  and  A.  Valentin,  J.  Mol.  Spectrosc. 

70,  18  (1978). 

R.  N.  Zare,  A.  L.  Schmeltekopf , W,  0.  Harrop  and  0.  L.  Albritton,  J.  Mol. 
Spectrosc.,  46,  37  (1973). 


T * 

• B « 

C ^ 

C H 

o ^ 


rs. 

CO 

CM 

CM 

ir> 

CO 

r>. 

VO 

a> 

o 

o 

VO 

JO 

rs. 

CVJ 

CO 

o> 

r~ 

CM 

CM 

CM 

Ch 

r~ 

VO 

in 

vD 

VO 

in 

CM 

p— 

in 

• 

* 

« 

• 

• 

• 

CO 

in 

CM 

CO 

o 

CM 

o 

in 

p— 

in 

VO  *— 

a%  00 

in  r— 

CO  o 

CO  o 

in 

CO  CO 

CO 

c*^  o 

CM 

VO  ^ 

^ CSJ 

CO  o 

in  VO 

00 

GO  O 
^ o 

^ o 

in  o 

VO  CO 

O)  o 

cr»  r- 

O r- 

in 

VO  O 

in  o 

r—  m 

r—  O 

in  o 

CM  in 

CO  o 

r-’  O 

vn  o 

CO  o 

CM  O 

f—  o 

00  •- 

CD 

CM 

00 

in 

p— 

in 

VO 

o 

o 

o 

r>* 

p» 

X 

X 

X 

CO 

CO 

00 

00 

CO 

CO 

CO 

CO 

> 

> 

> 

> 

> 

0) 

0) 

o 

CD 

o 

3: 

«c 

<0 

CM 

CO 

VO 

CM  ro 

a> 

o 

GO 

CM 

CM  *— 

vn 

O 

VO 

ja 

o 

CO  O 

o 

Cv 

OO 

VO 

in 

in  o 

VO 

a> 

CM 

VO 

in 

00  o 

VO 

in 

CM 

CO 

• • 

« 

• 

• 

« 

• 

CO  o 

in 

CM 

CO 

CM 

o 

CM 

CM 

o> 

pv 

CO 

VO 

vn  ^ 

o>  »— 

VO  O 

a»  o 

CM  VO 

in  tn 

VO  ^ 

oo  o 

o>  ^ 

^ r— 

O)  C) 

00 

^ CM 

o o 

00  CM 

O 1— 

in  ^ 

rv  CM 

vn  00 

00 

^ o 

VO  O 

CM  O 

vn  GO 

VO  O 

^ o 

GO  VO 

^ o 

CO  o 

VO  O 

in  o 

CO  f— 

^ o 

VO  O 

o o 

o 

ro  O 

r-  O 

if)  CD 

CM  O 

CM  O 

r-^  O 

rv  p- 

d o 

CM 

CM 

00 

p* 

CO 

VO 

p>— 

o 

o 

o 

pw 

X 

X 

X 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

> 

> 

> 

> 

> 

O 

a> 

0) 

C9 

GO 

a 

3: 

< 

< 

(d 

JO 

CO 

<Tt 

00 

p“ 

CO 

VC 

VO 

O) 

VO 

VO 

00 

00 

o> 

vn 

o> 

o» 

r- 

in 

CM 

ps 

VO 

o 

O) 

p» 

• 

• 

• 

• 

• 

« 

• 

vn 

CM 

CM 

o 

o 

CM 

O) 

p» 

VO 

o 

o 

o 

PM 

X 

X 

X 

o 

o 

o 

o 

o 

o 

o 

> 

> 

> 

> 

0) 

0) 

CD 

o 

< 

< 

Id 

XI 

assumed  computed  using  E 


